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Abstract The X-ray spectra of the magnet ar candi- 
dates are customarily fitted with an empirical, two com- 
ponent model: an absorbed blackbody and a power-law. 
However, the physical interpretation of these two spec- 
tral components is rarely discussed. It has been recently 
proposed that the presence of a hot plasma in the mag- 
netosphere of highly magnetized neutron stars might dis- 
tort, through efficient resonant cyclotron scattering, the 
thermal emission from the neutron star surface, result- 
ing in production of non-thermal spectra. Here we dis- 
cuss the Resonant Cyclotron Scattering (RCS) model, 
and present its XSPEC implementation, as well as pre- 
liminary results of its application to Anomalous X-ray 
Pulsars and Soft Gamma-ray Repeaters. 
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1 Introduction 

Anomalous X-ray Pulsars (AXPs) and Soft Gamma-ray 
Repeaters (SGRs) are a small class of slowly rotating 
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(5-12 s) neutron stars with emission properties much at 
variance with those of ordinary X-ray pulsars, both the 
young radio pulsars and the X-ray binary pulsars. They 
are called "anomalous" because their high X-ray lumi- 
nosity (10 34 — 10 36 erg/s) cannot be easily explained 
in terms of the conventional processes which apply to 
other classes of pulsars, i.e. accretion from a binary com- 
panion or injection of rotational energy in the pulsar 
wind/magnetosphcre. On the other hand, measurements 
of spin periods and period derivatives, when the lat- 
ter are interpreted as due to electromagnetic dipolar 
losses, suggest that these objects may host "magnetars" , 
i.e. neutron stars endowed with an ultra-strong mag- 
netic field (10 14 - 10 15 G, see Duncan & Thompson 1992; 
Thompson & Duncan 1993, 1996). The magnetar sce- 
nario appears so far very promising in explaining both 
the main energy source of these objects (the decay of 
the super-strong field) and the emission of the short, 
energetic bursts. Moreover, the magnetar scenario can 
account for the properties of Giant Flares, extremely en- 
ergetic transient events (L « 10 44 — 10 47 erg/s) detected 
from SGRs. However, alternative scenarios to explain the 
enigmatic properties of these sources have been invoked. 
Among these, models involving accretion from a fossil 
disk, left by the supernova event which gave birth to the 
neutron star, are still largely plausible (van Paradijs et 
al. 1995; Chatterjee, Hernquist & Narayan 2000). Very 
recently such a (maybe passive) disk has been indeed 
observed in the IR around AXP 4U 0142 + 61 (Wang, 
Chakrabarty & Kaplan 2006). Magnetars candidates are 
strong X-ray sources, and their X-ray spectra in the 0.5- 
10 keV band are usually fit with a two component model 
consisting of (besides interstellar absorption) a thermal 
blackbody with a typical temperature kT ^0.3-0.4 keV 
and a power-law with a relatively steep photon index, 
r ^2-4 (see Kaspi 2006 in this volume, and Woods & 
Thompson 2006 for recent reviews). In some cases, SGR 
spectra have been fit with a single power-law compo- 
nent, but recent longer observations showed that, also for 
these sources, a blackbody component is often required 
(Mereghetti et al. 2005). 
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Fig. 1 Modification of an initially Plankian spectrum 
(dashed line) by multiple cyclotron scattering for different 
values of to and (3 t h=0-3 (from Lyutikov & Gavriil 2006). 

Despite the fact that the blackbody plus power-law 
spectral model has been largely applied to magnetar 
spectra for many years, a reliable physical interpretation 
of these two components is still missing. Quite recently, 
an attempt to interpret the observed spectrum and long 
term variations of IE 1048.1-5937 in terms of a (non 
resonant) Comptonization model has been presented by 
Tiengo et al. (2005). Furthermore, the recent discovery 
of magnetar counterparts in the radio to the 7-ray bands 
(Camilo et al. 2006; Hulleman et al. 2000; Kuiper et al. 
2004) enforced the idea that their multiwavelength spec- 
tral energy distribution is by far more complicated than 
a simple blackbody plus power-law distribution (see also 
Den Hartog et al. contribution in this volume). 



2 The Resonant Cyclotron Scattering (RCS) 
model applied to magnetars spectra 

Following the original suggestion by Thompson, Lyu- 
tikov & Kulkarni (2002), it has been proposed by Lyu- 
tikov & Gavriil (2006) that the presence of a warm non- 
relativistic plasma in the magnetosphere of magnetars, 
might be responsible, through efficient resonant electron 
cyclotron scattering, of distorting the thermal X-ray emis- 
sion from the star surface. In particular, if a large vol- 
ume of the neutron star magnetosphere of the is partially 
filled by e ± -currents, the thermal (or quasi-thermal) cool- 
ing radiation emerging from the star surface will experi- 
ence repeated scatterings at the cyclotron resonance. 

The efficiency of the process is quantified by the scat- 
tering optical depth, r res . Following again Lyutikov & 
Gavriil (2006), and assuming a dipolar configuration for 
the magnetic field, this can be written as 

T~res = / O res ndl = Tq(1 + cos 2 a) 
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Fig. 2 The XMM-Newton spectrum of IE 1048-5937 fitted 
with an absorbed power-law plus a blackbody (top panel) and 
with the RCS model (second panel). 

where a res is the cross section for electron scattering in 
the magnetized regime, n is the electrons number density, 
a is the angle between the photon propagation direction 
and the local magnetic field, and 

7t e nr 

r ° = q ' 

6m e cu>B 

Here r is the radial coordinate, ujb = eB/m e c is the 
electron cyclotron frequency, and B is the local value of 
the magnetic field. At variance with the case of unmag- 
netized Thomson scattering, in presence of a magnetic 
field a res is resonant and is given by 

ot (1 + cos 2 a)uj 2 

where r = 4e 2 ui B /3m e c 3 is the natural width of the first 
cyclotron harmonic and <tt is the cross-section for 

unmagnetized Thomson scattering. The crucial point 
is that, in this model, the plasma distribution is spatially 
extended. Since the value of the resonant energy depends 
on the local value of the magnetic field, repeated scat- 
terings of photons could lead to the formation of a high 
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energy tail instead of a narrow line. At relatively large 
distances from the surface, the magnetic field drops to 
typical values such that scattering is resonant for soft 
X-ray photons, below ~ 10 keV. On the other hand, 
at these energies the resonant scattering optical depth 
greatly exceeds that for Thomson scattering, tt ~ max, 



Table 1 Best fit values of the spectral parameters obtained 
by fitting the XMM-Newton spectrum of IE 1048-5937 with 
the BB+PL and an RCS model (for more details about this 
observation see Mereghetti et al. 2004). Errors are at la con- 
fidence level, and the reported flux is absorbed and in the 
2-10 keV energy range. 
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where r e — e 2 /mc 2 is the classical electron radius. This 
implies that even a relatively small amount of plasma 
suspended in the magnetosphere of the neutron star may 
considerably modify the emergent spectra. There is now 
increasing support to the idea that the magnetar mag- 
netospheres might be twisted (Thompson, Lyutikov & 
Kulkarni 2002), in which case they may be partially 
filled by e~ currents with densities well in excess of the 
Goldreich- Julian density (which is expected in the case 
of a simple dipolar configuration, see Goldreich & Julian 
1968). 

A further effect is related to the fact that the spatial 
region, over which such currents are distributed, is per- 
meated by a highly inhomogeneous magnetic field. Al- 
though for soft X-ray photons the scattering conserves 
energy in the electron rest frame, thermal/bulk electron 
motion produces an energy transfer in the observer rest 
frame. Multiple resonant cyclotron scatterings of ther- 
mal photons will, on average, up-scatter the transmitted 
radiation, resulting in the formation of a non-thermal 
spectrum. The emerging spectrum will then appear dis- 
torted and, as it is the case for thermal non-magnetic 
comptonization, its shape may be well represented by 
the superposition of thermal and non-thermal compo- 
nents (see FigQ. 

Here we present an application of these models to 
magnetar spectra. Theoretical models are computed as 
in Lyutikov & Gavriil (2006), under the assumption that 
scattering occurs in a non-relativistic, static, warm medium 
and neglecting electron recoil, i.e. in the limit (3th < 1 
(where /3th, a model parameter, is the thermal velocity 
of electrons in units of the speed of light) and for low- 
energy photons, e <C mc 2 . Within the range of parameter 
values allowed by the numerical code, the non-thermal ef- 
fects are most prominent for models in which the plasma 
is mildly relativistic, (3 t h ~ 1, and the optical depth is 
reaching its boudary value r res ~ 10. 



3 The RCS XSPEC model 

In order to perform a quantitative comparison between 
fits to magnetar spectra performed with the model de- 
scribed above and with the canonical blackbody plus 
power-law model, we developed a code to upload the 
RCS model into XSPEC. We created a grid of intensity 
tables (through a Monte-Carlo simulation) for a set of 
values of the three model parameters, i.e. (3th, To and T, 
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0.03 

0.99 1.08 
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where T is the temperature of the seed thermal surface 
emission (assumed to be a blackbody). The parameters 
ranges used were 0.1 < (3 t h < 0.5 (step 0.1), 1 < r < 10 
(step 1) and 0.1 keV< T < 3 keV (step 0.05 keV). For 
each model, the spectrum has been computed in the en- 
ergy range 0.1-12keV (step 0.05 keV). The final XSPEC 
atable spectral model has therefore four parameters, the 
three listed above plus the last one being the normal- 
ization constant, which are simultaneously varied during 
the spectral fitting following the standard x 2 minimiza- 
tion technique. 

Note that a previous attempt to apply the Lyutikov 
& Gavriil (2006) model to spectral data was presented by 
the authors themselves. However, in that case (3th and To 
were fixed a priori during each fit and eventually varied 
manually. The fit itself was done only by varying T and 
the normalization factor, while in this case we are able 
to fit simultaneously all the four parameters, resulting in 
a more precise determination of their best fitting values 
and in a more reliable x 2 determination. 



4 Preliminary Results 

Here we present the preliminary results of the spectral 
modelling of magnetars with the Resonant Cyclotron 
Scattering model (RCS). We concentrate on the soft X- 
ray spectra (1-10 keV) of three sources, IE 1048-56, 
1RXS J1708-4009 and SGR 1806-20, and we use data 
obtained with the XMM-Newton satellite. Detailed in- 
formation about the observations, data reduction and 
analysis can be found in Mereghetti et al. (2004), Rea et 
al. (2005a) and Rea et al. (2005b), for the three sources, 
respectively. 

The fit of the X-ray spectrum of IE 1048-56 is pre- 
sented in Figure|2]and Table 1. In this case we find that 
the RCS model is successful in reproducing the data, as 
well as the canonical blackbody plus power-law model 
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Fig. 3 XMM-Newton spectra of 1RXS J1708-4009 (top 
panel) and SGR 1806-20 (bottom panel) fitted with he RCS 
model. The black-body plus power-law fits might be found in 
Rea et al. (2005a) and Mereghetti et al. (2005). 



(BB+PL). The two fits have the same number of de- 
grees of freedom. The value of the column density found 
with the RCS model is lower (although both consistent 
with Durant & vanKerkwijk 2006), while the tempera- 
ture of the thermal component T and the total flux are 
consistent with that found with the BB+PL modelling. 
These results are in agreement with the previous findings 
by Lyutikov & Gavriil (2006), who presented a similar 
attempt by using Chandra data and varying manually 
some of the model parameters to reproduce the spec- 
trum. 

As mentioned before, the technique used here is more 
reliable in providing a determination of the best fitting 
parameters. In this case, the source spectrum can be re- 
produced by invoking scattering by e~ with thermal ve- 
locity (3 t h — 0.41 and optical depth t = 1.9. 

On the other hand, the cases of 1RXS J1708-4009 and 
SGR 1806-20 are more problematic. These sources have 
much harder X-ray spectra (in the 1-10 keV band), that 
cannot be accounted for by the RCS model alone (see 
Fig.lSJ at least when the spectral deformation resulting 



from the resonant scattering process is computed in the 
simple way proposed by Lyutikov & Gavriil (2006). 

However, it is worth noting that, at variance with 
IE 1048-56, these two sources exhibit strong high-energy 
power-law tails which extend up to the 7-rays (Kuiper 
et al. 2006; Molkov et al. 2005; Mereghetti et al. 2005b). 
Motivated by that, we fitted the whole 1-50 keV spectra 
with the RCS model plus a power-law. For both sources, 
we obtained a good fit with a reduced % 2 of 1.08 and 1.1, 
respectively. Although extremely promising, these anal- 
ysis are very preliminary and they still require a careful 
assessment of intercalibration issues which arise when 
different detectors are used. Therefore, we do not dis- 
cuss the results in more detail nor we present here best 
fitting parameter values. These results will be presented 
in more details in a forthcoming paper (Rea et al. 2006, 
in prep.). 



5 Conclusion 

The preliminary results reported here show that the X- 
ray spectra of the AXP IE 1048-56 can be successfully 
reproduced by the RCS model. The quality of the fit 
is comparable to that obtained by using the canonical 
blackbody plus power-law model, with the advantage 
that the RCS model spectra provide a better physical un- 
derstanding of the presence of the thermal/ non-thermal 
components in the observed spectrum. On the other hand, 
this is not totally unexpected since this source has a 
non-thermal tail which is relatively soft with respect to 
other AXPs, and a similar result has already been found 
by Lyutikov & Gavriil, although their fitting technique 
was less accurate. For the other two sources considered 
here, 1RXS J1708-4009 and SGR 1806-20, we find that 
the RCS model alone cannot account for the relatively 
hard non-thermal component observed below ^10 keV. 

On the other hand, in both cases a composite model 
consisting of a RCS model plus a power-law can fit the 
X-ray spectra in the 1-100 keV band. This may suggest 
that, in the 6-10 keV energy range, the X-ray emission of 
these sources is contaminated by the hard non-thermal 
component which has beed recently observed with IN- 
TEGRAL at much higher energies (Kuiper et al. 2006). 
A similar hard tail has not been observed, at least so far, 
in the case of IE 1048-56. A detailed paper with a sys- 
tematic application of the RCS model to the whole class 
of magnetars and a study of their spectral and flux vari- 
ability in this scenario is in preparation, accompanied by 
a more complete discussion (Rea et al. 2006, in prep.). 
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